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DAMAGE STABILITY SYSTEM FOR THE SHIP HULL CHARACTERISTICS PROGRAM 


1. GENERAL 


The Damage Stability system for the Ship Hull Characteristics Program (SHCP) 
consists of two major supervisor programs; 


A. Subroutine CMPRTMNT 

This program supervises the generation of compartment offsets for com= 
partment bulkheads and all intermediate ship stations that lie within the 
bounds of the compartment, based on an outline description of each bulk- 
head. A Compartment Data Table (CDT), similar in form to the Ship Data 
Table (SDT), of offsets, sectional areas, centers of area, and breakpoint, 
is then created and stored on the Ship Data Tape (SDTPE) for later recall 
by Subroutine DAMAGE, 


B. Subroutine DAMAGE 

This program supervises the calculation of damage stability character 
istics of the ship minus selected flooded compartments for a range of draft, 
trim, and heel conditions. The sequence of events for each condition of 
compartments, draft, trim, and heel is as follows: 


1, The SDT is read from the SDTPE. 

2. The volumetric properties of the intact ship are calculated. 

3. For each selected compartment, the associated CDT is read from 
the SDTPE and the volumetric properties of the compartment are cal- 
culated and stored, 

4, The net volumetric properties and the righting arm for the ship 
minus compartments, i.e., the flooded condition, is computed and 


output for the condition is printed out. 


The Compartment Data Tables may be created by use of Subroutine CMPRT- 
MNT or by submitting compartment offsets in the same manner of submission as ship 


offsets for the Ship Data Table. If the compartments are submitted by the latter meth~ 
od, the LBP card must contain the Length oF the Ship, and the stations must be in terms 
of ship stations. 

C. Program Identification Numbers for calling by the Executive Routine 


CMPRIMNT....eeceeseesee O07 
DAMAGE.....-seeeccceees 008 


If. Subroutine CMPRTMNT 


A. 


Execution Time: 50 seconds/compartment, maximum 
Input: 
1. Ship Data Table (on tape) 


2. Card Types |, 2, 6: see NAVSEC documentation of the EXECutive Program 
for description of card types 1, 2, 6, E, which are always required. NAVSHIPS 
0900-006-5600. 


3. Card Type 25: Number of compartments; Format-15; NCOMP! 


For each compartment 

4. Card Type 26: Compartment serial number, compartment name, date, 
symmetry sentinal, number of bulkheads in the com- 
partment; Format-14, 10A4, 6X, 213; ISERNO, 
NAME 1S, DATE, |1SYM, NBULK. Symmetry sentinal 
definitions: 000, compartment on, or flooded on, both 
sides of the centerline; b-1, compartment on, or flooded 
on, the port side only; b+] compartment on, or flooded 
on, the starboard side only. (b represents blank) 


For each bulkhead 
5. Card Type 27: Number of outline points on the bulkhead, dis- 
tance of bulkhead from FP 9 (in feet); Format-I5, 
F1I0.2; NINT(I), XBLK(I) 


For each outline point 
6. Card Type 28: Z of outline point, Y of outline point; 
Format-2F 10.2; ZINT(I, IPT), YINT 
(!, IPT) 


Sample Input 


2 compartments, 2 bulkheads each, 3 outline points on Ist bulkhead of each 
compartment, 4 outline points on 2nd bulkhead of each compartment; 


Card deck by card type number: 1,2,6,25,26,27 ,28,28,28,27 ,28,28,28 ,28 
26,27 ,28,28,28,27,28,28,28,28,E 


SL ar 


e 


li! Subroutine DAMAGE 


A. Execution Time: (1+ number of compartments in set) x 5 x number of 
drafts x number of heels x number of trims. 


time is in seconds for a set of compartments 
B. Input: 
1. Ship Data Table and Compartment Data Tables (on SDTPE) 
2. Card Types!,2,6; see NAVSEC documentation of EXECUTIVE routine 
for description of card types 1,2,6 which are always required. 
NAVSHIPS 0900-006-5600 


3. Card Type 29: Number of sets of compartments (maximum :20); 
Format- 15; NSET 


For each compartment set 
4. Card Type 30: Number of compartments inset (maximum :25); 
Format- 14; LICOMP(!) 


5. Card Type 31: List of <erial numbers of compartment: in. set f 
and associated permeabilities (10 Compartments per card); 
Format= 10(14,F4.2); LACOMP(I,J), PERM8(I,J) 


6. Card Type 32: Number of trims, number of drafts, number of heels 
(maximum: 3,8,10} Format = 315; N4TRIM, N4DRFT, N4HEEL 


7. Card Type 33: Input trim angles in degrees (one per card); 
Format - F10.3; DITRIM 


8. Card Type 34: Input drafts in feet above baseline (one per card); 
Format - F10.3; DIDRFT 


9. Card Type 3%: Input heel angles in degrees (one per card); 
Format - F10.3; DIHEEL 
C. Sample Input 


Tw. sets of compartments, the first set with less than 10 compartments, 
the second set with more than 10, 5 drafts, no trims and 6 heels; 


Card deck by card type number: 1,2,6,29,30,31,31 ,30,31,32,34,34, 
34,34,35,35,35,35,35 ,E 


NOTE: There is no Card Type 33 because no trims were submitted. In this 
case a trim of zero is a’sumed and submission of Card Type 33 is 
not required. 


The Damage Stability Computer Programedescribed herein are for incorporation 
into the Ship Hull Characteristics Program developed by Mr. G. E. Hirschberg and 
Lt. R. L. Warters ot the Naval Ship Engineering Center, Washington, D. C. 


For documentation of those programs not part of the Damage Stability System, 
the reader should see the documentation produced ty NavSec: 


Nav Ships Number 


0900-06- TITLE 
5600 EXECUTIVE 
5510 AREAS 
5620 BA. 
5630 COEFIC 
5640 COEFI 
5650 DRAFT 
5670 FNTPIF, FNTP2F, FNTGIF, TRANSF, & TESTS 
5690 FI.NGTH 
5700 HYDRO 
5710 HILIM 
5720 OFFSET 
5730 STRNTH 
5740 TABLE 
5750 TESTI 
5760 TEST2 
5770 VOLUME 
5790 WAVE 
5800 WPLANE 
5810 WPLAN] 
5820 WETSRF 
5830 INTACT 


DOUBLE OFFSET MODIFICATION 


1. Modification Description: 


The modifications described herein permits YGHULL to handle all conceive- 
able hull forms with the exception of the trimaran for all of the major supervisors 
within the program system. 


The major modification to the program exists in Subroutine AREAS, Except 
for housekeeping changes such as changes in the DIMENSION and COMMON 
Statements and limits of some DO Loops the system remains unchanged. 


If the hull form is such that double offsets exist over some portion of the en- 
tire length = the catamaran hull form and a well are examples - the hull is sub- 
mitted to the program in two parts. One "hull" is the offsets of the outer portion 
of the hull. The second “hull" consists of offsets of the inboard portion of the 
hull. See Figure 1. The Subroutine AREAS then subtracts the properties of the 
"inner hull" from the "outer hull" to get the properties of the actual ship hull 
form. The properties of the section are then returned to Subroutine VOLUME 
and the system then proceeds normally. 


The basic operation of AREAS is now as follows: 


The program proceeds through AREAS normally to the point just before 
returning to the calling program which is Statement 1307. At that point the 
new variable ICAT is interrogated to see if it is equal to one, two or three. 
If ICAT is one, then AREAS assumes a normal ship and returns to the calling 
program. If ICAT is two, AREAS assumes a "dual hull" form and stores the 
properties computed so far in temporary accumulators. These are the proper- 
ties of the "outer hull". The variable iCAT is then set equal to three and 
the station indicator ISTAT is increased by NSTAT, the number of stations 
in one hull. The change in the index ISTAT is based on the premise that 
that two "hulls" are submitted in sequantial order with the "outer hull" 
being first, the “inner hull" being second and each hull having the same 
number of stations located af the same place along the longitudinal axis. 
The program then makes another pass through AREAS « + computes the properties 
of the "inner hull". 


This time when the program gets to the test of the variable ICAT following 
Statement 1307 the value of this variable is three (3). The program then 
branches to Statement 799 where the net properties of the section ave com- 
puted. ICAT is reset to two and ISTAT is reset to its original value. Control 
then returns to the calling progrom with properties of the net sections stored 
in the work table. 


It input: 


Each hull" is submitted in accordance with standard procedures, the last 
point on the last station of each hull having a breakpoint indicator value of 99999. 


SHIP wien @ WwWeELL 


| CATRMARAN 


eneelie 


OvTtR WELL INNER VOLE 


6 ElGvRe + 


DOUBLE OFFSET MODIFICATION (continued) | 


Each "hul!" must have data submitted at the same stations. If there fs no inner 
hull ata station - in the case of a well on inner huil exists only in way of the well 
- a data card must still be submitted for that station, the Y and Z values being zero 
with a breakpoint indicator value of 88888 or 99999 if it is the last station. 


On any given station, the data points on the inner hull may be independent of 
the data points on the outer hull. 


The outer hull is submitted first, the inner hull follows, and the LBP card is 
last. 


Card Type 1: The number 002 should be punched in cc. 51-53 if two hulls 
are being submitted; 000 for a regular hull. 


LONGITUDINAL BREAKPOINT MODIFICATION 


| Modification Description: 


The standard version of the Ships Hull Characteristics Program performs longi= 
tudinal integration from the first station to the last station by integrating groups of 
3 consecutive stations according to Simpson's first rule (1, 4, 1) and summing the 
results. If there are an even number of stations, or if the middle station in a group 
of 3 does not lie half way between the first and third, interpolation is performed to 
find the values of the missing, or evenly spaced, station. However, there is no 
method of taking into account any longitudinal discontinuities in the hull, such as 
the end of a poop or forecastle deck, or the boundaries of a well (the latter require 
incorporation of the Double Offset Modification to SHCP). The modification des- 
cribed herein eliminates this restriction. 


This is accomplished by submitting twice any station where longitudinal inte- 
gration should stop, and begin again. For example, in Figure 1 it ts desired to 
stop integration at the end of the sonar dome (and then begin again). Thus, Sta- 
tion 4 should be submitted twice. 


Figure 2 illustrates a slightly different case. Here the area of Station 5 is 
different when integrating aft from Station 1, than when integrating aft from Station 
5. Consequently, the first Station 5 submitted will describe the area 5A in Figure 2, 
and the second Station 5 submitted will describe the area 5B. The program integrates 
to the first Station 5 encountered (area 5A), and begins again at the second Station 
5 (area 5B), 


Figure 3 illustrates the use of this feature when using the Double Offset Mod= 
ification. Here a well begins at Station 6 and ends at Station 8. Thus it is desired 
to resubmit both of these station twice. They are submitted as illustrated 641, 6A2, 
6B1, 6B2, 8Al, 8A2, 8B1, 8B2. 


Each double station counts as two stations with respect to all arrays. This should 
be considered when counting total station so as not to exceed the maximum number 
permitted, which is 30. 


SUBROUTINE CMPRTMNT 


Description: 


This program supervises the generation of compartment offsets for compart- 
ment bulkheads and all intermediate ship stations that lie within the bounds of the 
compartment, based on an outline description of each bulkhead. 
The basic operation of the program for each compartment is as follows: 
1, The Ship Data Table (SDT) is read from the Ship Data Tape (SDTPE). 
2. The compartment title card is read. | 
3. The outline point ¢....,tion for all bulkheads are read. 


4. The outline point description of each bulkhead is tested to determine | 
if a double hull is required for the compartment. | 


5. Subroutine GENERATE is called which generates the compartment 
offsets. | 


6. The compartment offsets are transferred to the ship offset arrays. 


7. Subroutine OFFSET is called where the SDT for the compartment is 
created and stored on the SDTPE. 


BULKHEAD OUTLINE DEFINITION POINTS 


Each bulkhead is described by a series of Z and Y coordinates which de- 
fines the outline of the bulkhead section with reference to the ship section and 
the ships YZ coordinate system, These points shall be referred to hereafter 
as “outline points" or “OP's,” 


Outline points exist wherever there is a discontinuity in the section out- 
line, such as the intersection of a longitudinal bulkhead with a deck or the 
hull. Figure 1 illustrates the definition of outline points. 


Z 


Compartment 


Bulkhead X=Ship data table points. 


Y FIGURE 1 


It should be noted that the half section of the bulkhead is defined, the 
bulkhead being assumed to be a symetrical about the centerline. The value of 
the input variable 11SYM on the Compartment Title Card indicates which side 
of the centerline or both that the compartment lies on. 


The northeast most point of the bulkhead is OP #1. Successive OP's are 
taken at each discontinuity in the outline proceding in a clockwise direction. 
The last OP, #8 in Figure 1, is identical to OP #1, 

There are four basic types of OP's: 


1, Z=value, Y = value; describes an OP occurring entirely with the ship 
section. (OPs #1, 2,5, 6, 7 in Figure 1) 


2. Z=value, Y = 9999; describes an OP defining an intersection with the 
hull surface (OPs #3, 4 in Figure 2) 


3, 2= 9999, Y = value or 9999; describes an OP existing at the top of 
the ship section (Figuie 2, OPs 1, 3, 4) 


4, Z=~-9999, Y =value 9999, describes an OP existing at the bottom of 
the ship section. (Figure 2, OP #2) 


BULKHEAD CUTLINE DEFINITION POINTS (continued) 


a ¢ 
Op) 


Op#2 FIGURE 2 


Use of these four OP types must be strictly adhered to within the above 
definitions. Violations of these rules will produce errors in the bulkhead 
section definition. 


Figures 3 thru 6 further illustrate the use of the four OP types with a single 
hull vessel. 


The use of the four OP types is identical with a double hull vessel, Care 
should be taken that all OP values are defined with respect to the centerline 
of the ship not of the individual hull. tn other words, use the same coordinate 
system used to originally describe the ship. 


Figures 7 thru 11 illustrate principles of use with double hulls. 


Ore 
1 Value) Value; 
2 Value Value 
3 Value; 9999 
4 Value, 9999 
5 Values Values 
6 Value, 0 
7 Value7z 0 
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BULKHEAD OUTLINE DEFINITION POINTS (continued) | 
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FIGURE 6 
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BULKHEAD OUTLINE DEFINITION POINTS (continued) 


bd 
OP _s As 
1 Value] 9999 
2 -9999 Value 
3 Value3 Value3 
4 Value4-} 9999 
2 
BL 
CL 
Ship With Well 
FIGURE 7 
OP 2 ys 
a } 
1 Value, 9999 
2 Value 9999 
3 Value3 9999 
4 Value, 9999 
5 Values 9999 
BL 
| CL 
CATAMARAN 
FIGURE 8 
eo 2 3 
1 9999 9999 
2 -9999 9999 
3 9999 9999 
4 9999 9999 
BL 
cL 
FIGURE 9 “> 


NOTE: OP 2 is always required - i.e. lowest point on hull. 


BULKHEAD OUTLINE DEFINITION (continued) 


op gz y 45 
1 9999 Value, 
2 Values 9999 
3 Values 9999 
4 9999 Value, 
5 9999 Values; 2 
CL i 
FIGURE 10 
oP - ee 
1 9999 Value} 
2 Value 9999 
3 Value3 9999 
4 9999 Value, 
5 9999 Values—] 
CL BL 
FIGURE 11 


NOTE: The OP definition of Figures 10 and 11 are identical, even though 
in the former case the hull surface used is the "outer ship hull" and in the 
latter case it is the "inner ship hull." The program determines which ship 
hull to use by the following logic: OPs are taken clockwise from the NE 
point; therefore, if A Z between the first OP that is a hull intersection 
and the following OP which must also be a hull intersection, that is 23 - 
Zo, is negative, the ship hull must be the outer one; if 2 2 is positive, 
it must be the inner one. This follows from the definition of ship outer 
and inner hulls- the outer hull being that portion of the hull surfact lying 
outboard of the lowest point; the inner hull being that portion of the hull 
surfact lying inboard of the lowest point. 


COMPARTMENT D&FINITION 


Bulkhead B2 


COMPARTMENT DEFINITION (continued) 


Four "bulkheads" are required to describe the compartment illustrated on the 
previous page. Bulkheads 4 and C are obvious. At B two "bulkheads" are sub- 
mitted, Bl and 82, to alter the cross section of the compartment. These "bulkheads" 
have the same longitudinal location in the ship and are consequently treated as a 
double station. Thus, in integrating between A and B, the B) cross section is used; 
in integrating between Band C, the B2 cross section is used. 


The illustrated compartment may be described as one "compartment", i.e. from 
A-C, only when the entire run from A-C exists on both sides of the centerline, or 
on either side. If the symmetry of the compartment from A-B is not the same as that 
from A-C, i.e. A-B exists on both sides of the centerline but B-C exists only on one, 
the compartment must be submitted in 2 parts, or as 2 separate compartments. When 
executing the Damage Stability Program , it would then be necessary to call both com= 
partments, even though they are physically one. The sentinel 11SYM on the com- 
partment title card describes the compartment symmetry: 0 = both sides of centerline; 
+ = starboard side only; - = pe:r side only. The following table delineates a!l pos- 
sible combinations for the illustrated compartment. 


COMPARTMENT DEFINITION (continued) 
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SUBROUTINE CMPRTMNT 


START 


DO 900 JCOMP=1,NCOMP 
1F(JCOMP-1)20,20,10 


#1Q 


KK=1 


ZIS(1,J)=0. 
30 Aho 
Y (1,J)=0. 
M217 ")0. 


DO 13 J=1,10 


JBPIS(I,J) =0 
JBREAK(! ,J) = 0 


FLOWCHART 


NOTE: SDT: Ship Data Table 
SDTPE: Ship Data Tape 


Store ship serial number. 
Read number of compartments. 
Do for each compartment. 


IF JCOMP=1 ship data table 
has already been read in EXEC; 
GO TO 20. 


ISERNO=ship serial number. 


SDTPE exists (for OFFSET). 


Ship is on SDTPE (for OFFSET). 


Do not print out SDT (for OFFSET). 


Clear ship and compartment 
offset arrays to zero. 


Clear Ship and compartment 
breakpoint tables to zero. 


as 


SUBROUTINE CMPRTMNT 


ZINT(K,J) = 0. 
INT(K,J) = 0. 
"3 
& 


NP(I) = 0 
NPIS(I) = 0 

ng [XW=0. 
Coaun XIS(I) = 0. 


Call OFFSET 
O-" 
READ(60, 1000) SERNO , NAME1S, 
DATE, I TESTI ,INIT,11SYM,NBULK 
¥e@ NBA N 


READ(60, 1002) NINT(I) ,XBLK(I) 


READ(60, 1003)(ZINTI IPT), YINT(I,IPT), 
IPT=1, NPT 


#22 


FLOWCHART (continued) 


Clear compartment outline point 
arrays to zero, 


Clear ship and compartment number 
of points/ station tables to zero. 


Clear ship and compartment “station 
distance from FP" arrays to zero. 


Subroutine OFFSET read SDT of 
ship from SDTPE. 


Read compartment title card. 
Do for each bulkhead. 


Read number of outline points on 
bulkhead and distance of bulkhead 
from FP, 


Number of outline points on bulkhead. 


Read in outline points. 


SUBROUTINE CMPRTMNT 


DO 28 J=1,10 
#28 


BPIS(I ,J) = 0 
JCAT=1 


NPT=NINT(I}-1 


DO 24 J=2,NPT 


IF(ZINT(I,J)=ZINT(I ,J-1)) 
24 24,26 


IF(ABS(YINT(I,J))-.01) 
24 24,27 


#97 


24 


(Continue ) 


#400 
LSTAT = NSTAT 


° @ 


FLOWCHART (continued) 


Cleor compartment breakpoint 
table to zero. 


Assume compartment has single 
hull. 


Number of outline points on first 
bulkhead. 


Test to determine if compartment 
has double hull. Test each bulkhead. 


If the number of points on the 
bulkhead is equal to 3, the com- 
partment must have a double hull 
(See Notes #1). 


Otherwise, do not consider last 
outline point as it is the same as 
the first outline point. 


Test A Z between successive outline 
points. When A Z becomes posi- 
tive test the Y value of the outline 
point. If Y is zero the point lies 
on the ship centerline, which does 
not require a double hull; go on to 
next outline point. If Y is positive 
a double hull will be required; 
reset JCAT to 2. If any one bulk- 
head requires a double hull, so will 
the compartment as a whole. Go 
on to 400. 


Assume number of stations in com= 
partment equals the number of ship 
stations. 


SUBROUTINE CMPRTMNT 


Call GENERATE 


NSTAT=LSTAT 


ICAT = JCAT 
| LIMIT=NSTAT * ICAT 
DO 60 I=1,LIMIT 
| Nees) — 
X(EXIS() 
| iam DO 45 J=1,10 
é 


BREAK(! ,J}=0 


[K=1 | 


IDO 57 J=1,NPT 


| 


FLOWCHART (continued) 


Subroutine GENERATE creates 

the table of compartment offsets 

and stores them in the 1S arrays, 
The following statements transfer 
the compartment data to the ship 
tables. The compartment thus be- 
comes a "ship". Number of stations 
in "ship" equals the number of sta- 
tions in the compartment. 


Single/ double hull indicator. 


Array size for stations. 


Do for each station. 


Number of points/ station. 


Distance of station from FP. 


Clear ship breakpoint table to 
zero. 


K=index of compartment break- 
point table. 


L=index of ship offset table. 
J-index of compartment crfset table 
(do for each offset). 


Increase index of ship offset table 
by 1. 


24 


SUBROUTINE CMPRTMNT FLOWCHART (continued) 


If this is first offset go right to 
transfer area (454), 


If the index of the offset is less 
than the next breakpoint, go to 
transfer area (454). If the index 
is greater, ‘ncrease the breakpoint 
index (#56). 


IF(J-1-J BPIS(I,K-1)}54,53 ,54 
0 


If the index is equal to that of the 
next breakpoint, test to see if the 
previous point was also a breakpoint 


#53 If it was not, go to transfer area 
(#54). If it was, two successive 
breakpoints have occurred, and a 
space in the table must be left for 
insertion of an intermediate point. 
L=i+1 Insert the signal value of Z=999999 .0 
and increase the "ship" offset table, 
#59 L, by 1. 
Y(I,L}YIS(I,J) 
Z(1,L)=ZiS(I,J Transfer the offset in the compart- 
ment table to the ship table. 
1F(J-JBPIS(I ,K))57 ,55 ,57 ae If this offset was a breakpoint, store 
' the index in the "ship" breakpoint 
5 table. Go on to next offset 


#56 
K = K+] 


Increase index of compartment 
breakpoint table by 1. 


Number of points on station = final 
value of L. 


25 


wantenmnee as new xn amet 


SUBROUTINE CMPRTMNT FLOWCHART (continued) 


Station number of compartment 
stations expressed as a % of LBP. 


STATNO(I) = X(1)/ XLBP 


Put "ship" on tape (for OFFSET). 
For OFFSET. 

Print out "ship" data table (for 
OFFSET). 


| 
Create SDT for compartment and | 
write on SDTPE, also print SOT. | 

| 


Restore ship serial number. 


Return to EXEC(YGHULL). 


26 


SUBROUTINE GENERATE 


Program Description: 


GENERATE supervises the formulation of offset tables for the compartment 
based upon the general bulkhead descriptions and the ships offsets as delineated 
in the Ship Data Table (SDT). 


The operation is basically a merge of bulkheads and intermediate ship 
stations contained in the SDT into a table giving a full description of the com- 
partment. This manipulation involves three basic sets of arrays: 


1. Bulkhead arrays. 


XBLK( ) The distance of each bulkhead from the F.P. 
ZINT( , ) The Z value of the bulkhead outline points. 
YINT( , ) The Y value of the bulkhead outline points. 

2. Ship Data Table 


X(_) The distance of ship stations from the F.P. 
Z(_, ) The Z coordinates of the station offsets. 
Y(_, ) The Y coordinates of the station offsets. 


3. Compartment Data Trble 


XIS(_) The distance of the compartment station from the F.P. 
ZIS( , ) The Z coordinates of the station offsets. 
YIS( , ) TheYcoordinates of the station offsets. 


GENERATE works only with the X _—_ arrays, utilizing subroutines SECTN, 
STORE, and MOD for modifying the Y__and Z__ arrays. In essence,it determines 
the station indices in the 1S arrays where the bulkheads and intermediate ship 
stations will be located. Figure 1 illustrates the transfer of information involved. 


— 


SUBROUTINE GENERATE 


KSTAT 
JSTAT 


LBULK 
NSTAT 
LSTAT 

MSTAT 


ICAT 
JCAT 


eels 


SYMBOL LIST 


«> 


Inciex of ship stations in the Ship Data Table. 


Index of compartment stations in the Compartment 
Data Table. 


Index of bulkheads in bulkhead arrays. 
Number of stations in one ship hull. 
Number of stations in one compartment hull. 


Index of Ist station aft in X, Y, Z array of foreward 
most bulkhead. 


1 = ship single hull; 2 - ship double hull. 


1 = compartment single hull; 2 - compartment double 
hull . 


SUBROUTINE GENERATE FLOWCHART 


Set bulkhead index LBULK, 
compartment section index 
JSTAT, and control character 


IFCT equal to 1. 
K— Do 944 ISTAT = 1,NSTAT 
| IF(X(ISTAT)-XBLK(1))}944,944,9010 Locate first ship station aft 
| of foreward most bulkhead. 
KSTAT = ISTAT Store index of this station 
MSTAT = KSTAT in KSTAT and MSTAT. 
D 
to 7 
yes 
A | F(LBULK-NBULK) If this is the last bulkheed, 
20,20, 100 no go to #100, 
So 
INBLK(LBULK) = Store index of bulkhead in 
JSTAT-LSTAT*(IFCT=1) compartment section arrays 


in INBLK (LBULK). 


Store X of bulkhead in XiS 
array. 


JSTAT = JSTAT + 1 


Increase index of IS arrays 
by |. 


SUBROUTINE GENERATE FLOWCHART (continued) 


If this is last bulkhead, 
go to 100 


If the ship station at X(KSTAT) 
is foreward of the next bulkhead 

XBLK (LBULK+1), this ship station 
lies within the bounds of the 
compartment. Transfer the off- 
sets of the ship station to the 
sections position in the com- 
partment (IS) arrays using Sub- 
routine STORE, store the station 
location in XIS(JSTAT); and 
increase the compartment section 
index, KSTAT. 


If X(KSTAT) is aft of the next 
bulkhead X(LBULK+1), increase 
the bulkhead index by 1, go to 


LBULK = LBULK+1 
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SUBROUTINE GENERATE FLOWCHART (continued) 


#100 
GO TO(110,110,220) ,ICAT If this is the 2nd pass through 
F 1,2 for a double hull ship, GO TO 220, 
110 
LSTAT=JSTAT=1 LSTAT = number of stations, 
including bulkheads, in the 
compartment. 


IF(LSTAT.GT.2) GO TO 150 
If LSTAT is not greater than 2, 
no station in the Ship Data 
Table lay between the two bulk- 
heads. An intermediate station 
must be generated in order to 
preserve the curvature of the hull. 


LSTAT =3 Number of compartment station 
=3, 
XIS(LSTAT)=X15S(2) Move 2nd bulkhead to station 3. 


XI1S(2) = (XIS(LSTAT) + XIS(1))/2. New intermediate station halfway 
between bulkheads. 
INBLK(2)=3 2nd bulkhead located at station 3. 


JSTAT =2 Station index = 2. 


CALL SECTN Compute offsets of intermediate 


station in Subroutine SECTN. 


JSTAT = 4 Station index = 4, 


) 


SUBROUTINE GENCKATE FLOWCHART (continued) ~ 
. 
FCT =3 Control character. 
f{ ~ ~ © 00 1751 =1, LSTAT 
| 
Store X of compartment stations 
| LL=I+LSTAT in positions where inner com= 
partment hull will be stored. 
| t 
175 Storage is sequential = outer 
lg — — -{XIS(LL) = XIS(I) hull, then inner hull. 
|F(IFCT.EQ.3) GO TO 230 
i GO TO (230,210,220), ICAT re) If ICAT = 1, the ship has a | 
single hull, GO TO 230, | 


2 
If ICAT = 3, the ship has a 
double hull and this is the 
£210 . second pass, GO T0220. 


If I\CAT = 2, the ship has a 
double hull and this is the 
first pass. The inner hull must | 
now be considered. Set 

LBULK=1 | ICAT = 3 and IFCT = 2 for 
contro! purposes; reintialize 
LBULK, the bulkhead index, 
and set the ship station Index 

KSTAT at the location of the 
first Inner hull station aft of 
the foreward most bulkhead. 

. Go back to store Inner hull 
KSTAT = MSTAT+NSTAT ship stations, 


SUBROUTINE GENERATE 


#235 
KSTAT = 1 
JSTAT= 1 
Call STORE 


Call MODIFY(IMOD) 


FLOWCHART (continued) 


Reset ICAT 


_ Indicator for Subroutine MODIFY 


signifying that bulkheads, as 
distinguished from intermediate 
stations, are being considered. 


Index of Ist bulkhead to be 
calculated by Subroutine SECTN 


If the first bulkhead lies on 

the first ship station, Subroutine 
SECTN may not be used. 

The offsets of the first bulkhead, 
which Is the first station, are 
transferred to the IS arrays via 
Subroutine STORE. Subroutine 
MODIFY is then called to modify 
these offsets in accordance with the 
bulkhead outline description con- 
tained in the arrays ZINT and YINT. 


Low is reset to 2, the second 
bulkhead. 


SUBROUTINE GENERATE 


Do 245 LBULK = LOW, NBULK 


JSTAT = INBLK(LBULK) 


Call SECTN 
all MODIFY (IMOD) 


FLOWCHART (continued) 


Do for each bulkhead. 


INBLK(LBULK) contains 
the index in the IS arrays 
where bulkhead LBULK is 
to be stored. 


Subroutine SECTN calculates 
the offsets at the bulkhead 
location of the full ship 
section and stores them in 
the !S arrays at location 
JSTAT 


Subroutine MODIFY modifies 
these offsets in accordance 
with the bu!khead outline 
description contained in the 
arrays ZINT and YINT. 


SUBROUTINE GENERATE FLOWCHART (continued) 


LBULK = | Bulkhead index is set to |. 


Indicator for Subroutine MODIFY 
indicating that intermediate stations, 
as distinguished from bulkheads, are 
being considered. 


Do 270 JSTAT = 2,LSTAT Go through all compartment stations. 


If JSTAT is less than |NBLK(LBULK+1), 
which fs the index of the next bulkhead, 
then the station JSTAT in the IS array 
is an intermediate station; Subroutine 
MODIFY modifies these offsets in ac- 
cordance with the bulkhead outline 
description contained in the arrays 
ZINT and YINT. If JSTAT = INBLK 
(LBULK+1), the station is a bulkhead 
and has already been operated on (See 
previous page). If JSTAT > INBLK 
(LBULK+1), increase LBULK by 1. 


1 F(JSTAT-INBLK(LBULK+1)) 
260,270,255 


Call OUTPT Print offsets of compartment. 


Return to Subroutine CMPRTMNT 


wwe oe Tae 


SUBROUTINE SECTN 


Description: 


Generally, hull offsets are given only for ordinate stations. Many appli- 
cations, however require offsets at intermediate locations. SECTN finds the 
hull offsets for an intermediate section lying between two ordinate stations. 


Fig. 1 shows a typical situation with an intermediate section to be defined 


wW.T. | 

"7" BULKHEAD | 
| 

| 


PanENY Ky Ke K, $ STATION 
FIG. 1 


between ordinate stations Ky and K3. Since the intermediate section (abbr. 
1,5.) is taken closer to ordinate Ko than K3, it is reasonable to assume that 
the hull at |.S. will be shaped more like the section at station K> than that 
at Kj. Therefore if offsets for Ko are given at certain drafts Z(Ko,1), l=1, 
NP and we assume this set of offsets properly describes the hull shape, then 
defining offsets for |.S. at these some drafts should give a reasonably good 
representation of the hull shape at !.S. This means that an, irregularities 
in Ka will show up also in !.S. To further emphasize the similarity, if sta- 
tion Ka has a breakpoint at draft Z(Kp ,1BP) we define a breakpoint for |.5. 
at this draft. 


Having chosen a draft, Z(Kg,1), offsets on K, and K3 can be found either 
because they themselves are given as Input or can be interpolated from the 
offsets which are. When possible, then, a smooth second order curve Is passed 
through the points on K;, Kp and Kg defined by the given draft and the corres- 
poinding offsets. The offset on 1.S. for a draft of Z(Ko,1) is then assumed to 
Ite on this curve, its exact value being determined from the relative spacing 
(measured In distance from the F.P.) of K), Ko, 1.S., K3. 


SUBROUTINE SECTN (continued) 


The meaning of the "when possible is illustrated in Figs. 2., 3. 


4 t Ks s 4 ji STATION # 
? Ky Ks 


NsTaT ETAT NETHER a 4 sTaTion % 


Figure 2 illustrates the situation when K, and Ko both designate the same 
"double" station which fs the usual manner in which abrupt changes in the hull 
(discontinuities in the waterplane or half-breadth plan) are handled. Figure 3 
illustrated an intermediate section located adjacent to the first and last stations. 
In all three cases and in any of the possible combinations there are not three sta- 
tlons which can give a meaningfully continuous curve with which to interpolate 
the |,S. offsets. In these cases a simple linear interpolation is used. 


When all the required offsets have been stored the intermediate hull section 
is defined unless the ship has a "double hull". In this case only the outer hull 
has been defined and a second pass through the routine is necessary to handle 
the inner hull. 


« 


37 


SUBROUTINE SECTN (continued) 


Program Method (Introduction) 


Below is a somewhat simplified functional diagram which should help to 
clarify some of the indexes and operations used in SECTN. 


Ineyr 


x xfs 

kosrarns| | mecrstarie | 

[nsrats | tf rstate | i 

NP NI2S 

= Ee ao 
Ks | Ka] rereme | | 

Y YLS 


= EE 


OvurPeut 


es ee 
Z constructed like Y table. ZIS table constructed like YIS table. 
JBREAK constructed like Y table JSBPIS table constructed like YIS table 
except bottom row index = 10, except only has 10 rows, the maximum 
the maximum number of break- number of breakpoints, instead of NPIS 


points, instead of NPT. rows. 


baeetmamatne se eres sere ee ee 


SUBROUTINE SECTN (continued) 


1) Basically, a single pass through SECTN selects 3 indices K) Ka Kg (the 
notation is somewhat altered in the program listing) representing consecutive 
ordinate stations such *hat !.S, is closest to ordinate station Ky. These three 
indices specify three numbers in the NP array representing the number of offseis 
on each station. The indices also specify three columns in each of the Y and 

Z arrays. SECTN operates on these columns in such a way as to create a new 
column 


NPIS | NPIS 
and a new column ZIS(JSTAT,JJ) | in the YIS and ZIS 
JJ=1 i JJ=1 


YIS(JSTAT, JJ) 


arrays respectively. The choosing of which column is to be filled in, i.e. the 
setting of JSTAT, is external to the routine and need not concern us here. The 
entries in these columns represent the offset coordinates for !.S. Also, the num- 
ber of offsets NPIS(JSTAT) is entered from 


10 
Nata] and the breakpoint table JBREAK(KSTAT,JBRK 
STAT=Ko JBRK=1 


0 KSTAT=Ky 
is transferred to JBPIS(JSTAT ,JBRK Thus all the data concerning |.5. 
JBRK=1 


is entered in the IS arrays. 


2) This describes a basic run through the program. The only time this descrip- 
tion is not sufficient is in the case of a double hulled ship. In this case a second 
pass through SECTN to create the entries for the inner hull is necessary. For a 
double hulled ship the Y, Z, NP, and JBREAK arrays are each expanded by a 
factor of 2 to 2*NSTAT in the horizontal dimension. This is because at each 
station there are now twice as many entries (inner as weil as outer hulls). All 
the entries for the outer hull precede all the entries for the inner hull. 


Similarly, in SECTN, all the IS arrays (except XIS since the locations 
of the sections remain the same for inner and outer hulls) have their horizontal 
dimensions doubled to 2*LSTAT. 


SUBROUTINE SECTN (continued) | 


On this second pass through SECTN, the proper data for station K2, for 
example is now found NSTAT columns to the right in column Ky + NSTAT. 
Likewise, if the |.S. offsets on the outer hull were stored in stains NSTAT 
of YIS and ZIS the offsets on the inner hull will be stored at columns JSTAT 
+ LSTAT. 


eed 


Now follows a detailed description of the actual! progrom. in reading 
this a listing of SECTN should be referred to, 


Detail: 


(Variable names are defined in the next section) The upper limit, LIM, 
of a DO statement is set to the number of ordinate stations !ess one. The func- 
tion of the loop is to locate the position of 1.S. between two adjacent stations. 
To do this the following function is calculated: 

DIS=(X;4.-XI1S) 

Bak sa | 


The properties of DIS are given in the table below for XIS#X(1). 


TABLE | 


DIS= Relation between arguments 


<0 XIS<Xj<X;+1 or XX, X18 

=0 XIS = Xj+] 

70 XjeXI*X iy | 

>) —— (K1S=X;)<(Xj;4.1=-XIS) i.e. XIS is closer to X; 
<| (Xj4+1-X1S)<(XIS=X;) je. XIS is closer to X;4} 


e2a0oao a 


The loop is repeated until either condition b. or c. occur. In the event 
DiIS=0 the |.S. lies on one of the ordinate stations, namely on station JW+1. 
The program exits from the loop and calls a subroutine. Subroutine STORE, DIS=0 or 
which transfers entries from the NP, Y, Zand JBREAK arrays into the corres=- condition b. 
ponding |S arrays as described under 1) above. In this connection it may be 
useful to think of the index JW+1 as corresponding to the index Ky. 


ET ae a pee 


SUBROUTINE SECTN (continued) 


After returning from STORE the program has finished unless the ship is 
double hulled. If so ICAT=2* and a computed GO TO sends control to State- 
mett 230. Here ICAT is set to 3 to prevent a third pass through SECTN and 
JW, which for this case corresponds to Ko-1, and JSTAT is transformed as in 
2) above. Control is transfered to Statement 5 which, since DIS=0, again 
results in a CALL STORE and then a transfer to Statement 240 (because ICAT 
=3). Statement 240 resets JSTAT and ICAT for use in the calling program and 
returns control to the calling program. 


In the event DIS 0, 1.S. lies between ordinate stations JW and JW+1. DIS > 0 oF 
The program exits from the loop and checks for conditions d. ore. of Table!. condition c. 


NOTE: From here on the formula 'A¢—>B' means that the symbol, A, occurring 
in the program listing corresponds to the symbol or phase, B, occurring in the 
se-tion of this documentation called "Program Description". 


The object at this point in the program is to designate a variable NCLOS 
Kg , and a variable NFAR Kg (see Fig. 1) unless 1.5. lies to the right of 
Kg in which case NFAR K). In the first case X(NFAR) X(NCLOS), in 
the second case the inequality is reversed (Statements 30 to 40+00001). In 
either case X(NCLOS) is the !ocation of the closest ordinate station to 1.5. 
Statement 50 now interpolates, in a straight line sense, the height of the keel 
of |.S., using statement furiction FLINT. i.e. The lowest point of |.S. is 
considered to lie on the line joining the lowest points of stations NFAR and 
NCLOS 


We know at this point between which two stations the 1.S. lies and, through 
the indices NCLOS and NFAR, which of the two is closer. 


A) The Ideal situation at this point from the point of view of accuracy 
is to have the |.S. positioned as shown in Fig, 4 (identical to Fig. 1) or 5. 


! 
eae ~ 
NBSTAT = ATAT-4 NEAR NCLos WCLOB- 4 a ‘ STATION # 
ey) (na) (aD 


*ICAT ts an indicator set by the calling program. 


SUBROUTINE SECTN (continued) 


FIGS 


wsvar NSTAT-4 a ' STATION & 


In these cases, to find the coordinates (ZIS(JSTAT,1), YIS(JSTAT,1)) of 
point Q, which are after all the quantities we are looking for, the procedure 
is as follows (see loop DO 150, Statement 130+00004): 


ZIS(JSTAT,l) is simply the draft, Z(NCLOS,|) of Pi. YIS(JSTAT,1!) must 
be interpolated from the offsets of points p!, Pi, and Pi. The drafts at 
these points are all the same. The offsets at each station are found by 
calling WPLAN] which returns these offsets in the array Yl. WPLAN1 
is thus called 3 times per point Q, once for each station. (DO 140, 
Statement 130+00006). Then the second order curve, whose coefficients 
have been found in Statement 130°(CALL CQEF11) based on the ordinate 
values of the sections, is "passed through" Pi, Pi, Pi. (Statement 150). 
This process is repeated for each point on station number NCLOS (i.e. 
N2). 


This discussion distorts slightly one minor point. It will be noticed from 
the listing that the offsets and drafts for the ordinate stations have a different 
index, JL, from the |.S. which uses JJ. This differs from the above usage which 
seems to imply a one-to-one correspondence between points on station NCLOS 
and points on 1.S. This would indeed be the case except that, whereas points 
in the output array YIS always start at the beginning with the index JJ=1, it 
is possible that there may be some unusable points on the ordinate station NCLOS. 
This will always happen if, for example, station NCLOS has a point whose draft 
is less than the lowest draft on station NFAR. There would than be no way of 
interpolating an offset on NFAR as explained in the preceding paragraph. 
Therefore, the indexing of the ordinate station points begins with the first usable 
point, JIST, on station NCLOS (lower index, DO 150, Statement 130+00004). 
JIST for the cases pictured In Figs. 4 & 5 is calculated, in Statement 110 through 
120, as the first point of closest station with a draft greater than the drafts of 
the first points of either neighboring station (M1 or M2). 


*NOTE: The statement: 130+00001, 2, 3 were inserted to prevent the coefficients 
returned from COEFII from being written over. The coefficients are used in 
Statement 150, but COEFII is again called in Statement 140+00001 by Subroutine 
WPLAN1. 
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SUBROUTINE SECTN (continued) 


B) The preceding discussion holds almost unchanged for the cases pic- 
tured in Figs. 6 and 7 in which it has not been possible to use a station 
to the right or left of NCLOS (resp. for Figs. 6 and 7). 


FIG © 


a ' STATION K 


(eas) 
This situation arises most frequently when NCLOS is a double station, but 

can also occur when NCLOS is either the first or last station on the ship. 

(It is implicitly assumed here that no great error is introduced in the inter- 
polation for Qdue to 1.S. being near one end or the other of the fitted curve). 


C) The only remaining cases are those pictured in Figs. 2 and 3 where it 
is not possible to use any ordinate station to the right or left of NCLOS and 
NEAR and a linear fit is necessary. 


For these cases Statements 160 to 170 perform an analogous function to 
Statements 110-120 discussed under A) above. Similarly, 180 to 200 are com- 
pletely analogous to Statements 130+00004 through 150. The only difference 
is in the number of stations and the formuia used for the interpolation. 


Once the offsets for |1.S. have been calculated, control is transferred to 
Statement 220 which is a DO which transfers the breakpoint table, JBREAK for 
NCLOS into the JBPIS array. Following this NPIS(JSTAT) = the number of data 
points on 1.5. is set equal to NP(NCLOS), the number of data points on NCLOS. 


At this point all the necessary arrays have been filled and the program 
checks for a double hull! (cf. para. 2), P.4). 
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SUBROUTINE SECTN (continued) 


The only statements remaining which have not been covered in detail are 
50+00001 through 100 which set the indices NI, N2, N3, M1, M2. But | think 
that the logic and the output of this section are sufficiently clear from the figures 
to warrant only this remark: that a statement of the form I F(X(L)-X(L+1))A,B,A 


where L=NCLOS or NFAR, is a check for a double station as discussed in paragraph 
B)P. above. 


SUBROUTINE SECTN SYMBOL LIST 
| 

DIS Function used for checking location of |.S. 
JIST Index of first usable point on NCLOS for interpola- 

tion. 
JW Index of lower ordinate station of interval containing 

| e S e 
LIM Number of ordinate stations less than 1. 
N1,N2,N3 Indices for curve fit, N1 corresponds to rightmost 


station used, N2 the middle station, N3 the leftmost. 
(F.P. is at right of profile). 


NCLOS Index of ordinate station closest to 1.S. which fs al- 
ways one of the stations used in the curve fit. 


NFAR Index of ordinate station next closest to 1.8. which 
is also always used in the curve fit. 


Mi, M2 Corresponds to two stations used in curve fit other 
than NCLOS. M1 is always less than M2, 


NOTE: ‘Variables, apart from general indices, which appear in the program listing 
and do :.ct uppear above are treated elsewhere in the YGHULL documentation. 


SUBROUTINE SECTN MACRO FLOWCHART 


Find interval containing 


compartment station 
JSTAT (5) 
Does the 
Ys closest station lie 
on the bulkhead? 
Ne 
(Ss 20 Sie png ary heels TY nl eat = | 


The other station of 

the interval containing: 

JSTAT is indeed with | 
NFAR 


Find closest station 
to compartment sta- 
tion JSTAT 

----» NCLOS 


tore ship station 
JW+1 parameters 
in arrays for com= 
partment station 


60 | 


interpolate keel draft and 


eas) halfbreadth for C.S. 


Is closest station 
Foawap aft or f- dof C.S. 


\ ISim1? 
85 a Jo 
Is the closest 


AeTt 


Is closest ship 
statior the first station [Ne 
: ship? 


double station? 


Is station 
near a double /& 
station? 


Is station 
near a double = 
station? 


SUBROUTINE SECTN (continued) MACRO FLOWCHART 


@) & O & 


Entrance through either of these 


statements (90,100) means that a Since the subroutines which 
linear interpolation {s required, perform the hal febreadth 

as tracing the preceding logic interpolation expect to have 
{will show, the stations given from fore to 


aft, the indices NCLOSE, NFAR 
must be arranged tn ascending 


160 numerical order, in addition, a 
Find first usable tiiird station is selected and its 
point on station Index stored in numerical order. 
NCLOS |} EO epee > NI,N2,N3 

e--> JST . The indices of the two neighbor 

stations NCLOS and NFAR are 

180 stored also in numerical order 
Perform linear | = [| ssese > M1,M2 


interpolation for 
ompartment sta- 
ion offsets. 


(30 


Perform quadratic 
interpolation 


190*0000)/, 200 
19090000, YT.) 
210 


Store 
parometers. 


Modify 
Indices. 
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SUBROUTINE STORE 


Purpose: 


STORE fills the compartment section arrays (__1S arrays) when the section corres- 
ponds to a station in the Ship Data Table. 


Description: Input consists of: 


a. The number of points on the station to be transferred. 
b. The offset tables (1.e., Y and Z arrays) from the Ship Data Table. 
c. The breakpoint table, JBREAK, also from the SOT. 


Output consists of the same numbers transferred to the 1S arrays NPIS, YIS, 
ZIS, JBPIS (see diag. 1 of SECTN documentation). 


The transfer is accomplished by two arithmetic statements and two DO LOOPS. 
Those loops operate on the offset tables and breakpoint table in that order. 


SUBROUTINE MODIFY 


Description: 


This program takes the ship section and modifies it in accordance with 
the description given by the bulkhead outline points. For example, Figure | tllus- 
trates the ship section and the associated table of offsets; Figure 2 illustrates the 
bulkhead outline and the associated table of outline points, Figure 3 illustrates the 
resulting compartment section description and the associated table of offsets. 


FIGURE 1 


x 


| indicates outline point 
Y 


FIGURE 2 


a See ees e 


SUBROUTINE MODIFY (continued) 


OUTER HULL 


| Indicates offset 


x 


FIGURE 3 | 


re 


SUBROUTINE MODIFY SYMBOL LIST 
| Arrays 


1, Compartment Data Table arrays 


XIS() La distance of the compartment station from the 
YIS(, ) The Y coordinates of the station offsets. 

ZIS( , ) The Z coordira tes of the station offsets. 

NPIS( ) The number of offsets on each station. 

JBPIS( , ) Table of breakpoints for each station. 


2. Temporary storage arrays for ship section. 


YIST( , ) The Y coordinates of the station offsets. 
ZIST( , ) The Z coordinates of the station offsets. 
NPIST( ) The number of offsets on the station. 
JBPST( , ) Table of breakpoints on each station. 


3, Bulkhead Outline Point arrays 


XBLK( ) The distance of each bulkhead from the FP. 
ZINT( , ) The Z value of the bulkhead outline points. 
YINT( , ) The Y value of the bulkhead outline points. 
NINT( ) The number of outline points on each bulkhead. 


I! Indices and Control Sentinels 


INDX1 Outline point index. 

INDX2 index of outline point with lowest Z. 

INDX3 Danke Hull Index for _ IST arrays; l=outer hull, 2=inner 
hull. 


SUBROUTINE MODIFY 


INDX4 
INDX5 


INT 


INT] 
JSTAT 
LMT 


LMT1 


LBULK 
L2 
L3 
L4 
L5 


Lé 


SYMBOL LIST (continued) 


Index of outline point just below (lower Z) INDX!. 
Outline point index. 


Lower limit of DO LOOP for transferring ship section 
offsets to Compartment Data Table. 


Lower limit of DO LOOP for outline points. 
Station Index in Compartment Data Table arrays. 
Upper limit of DO LOOP for outline point. 


Upper limit of DO LOOP for transferring ship section 
offsets to Compartment Data Table. 


Bulkhead index, 

Offset index in Compartment Data Table arroys. 
Index of intial ship section breakpoint to be tested. 
Breakpoint index in Compartment Data Table arrays. 


Control sentinel: if outiine point being considered 
Is lowest point, 2 if ary other outline point. 


Control sentinel: 2 if outer hull, 1 if inner hull 


SUBROUTINE MODIFY (continued) 


NOTES REFERENCED BY FLOWCHART | 


Einal Compartment Section 
b 

work (only ZIS | 
up for shown) | 
outer 
hull 

outer 

hull 

Inner 
work hull 
down for inner hull 

__IST arrays __INT arrays __IS arrays 
1. The DO 500 LOOP goes from outline point to outline point (INT) merg- 


ing with the ship section offsets (_IST) to form the Compartment Section Offset Table 
(_1S). This table Is created starting at the lowest point, consequently, the Ship Sec= 
tlon Offset Table must be looked at from the lowest point and the Outline Point Table 
must be looked at starting at the lowest point. However, this latter table is in order 
from the Northeast most point (See Bulkhead Outline Definition Points), clockwise, 
and consequently the lowest outline point lies somewhere in the middle of the table. 
Thus the lowest outline point is located and its index stored in INDX2. Since the 
DO LOOP Index must run down the table, the index of the outline point used in ex- 
pressions must be altered to increase in the order Increasing Z's. 


When working on the outer hull (using the above illustration) the DO LOOP 
Index runs from 1} to 4 while the actual Index should run from 4 to 1. The Index ex- 


pression becomes 
INDX1 = LMT =! +1 
where LMT Is the number of points In the table (4) and | is the DO LOOP index. 


While working on the Inner hull the actual Index corresponds to the DO 


LOOP Index 
INDX) =I. 


Rather than use a COMPUTED GO TO to determine which index ts valid, 
a single expression is used. 


* eieiaaneiieniaise. ~* en i nce AC PS STON = EG SE TI Gs IE AAAI CIC, SOT = 
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SUBROUTINE MODIFY (continued) 


INDX1 = (LMT-I+1) * (L6-1) +1 * (2-L6) 
where L6 is 2 when working on the outer hull, and L6 is 1 when working on the inner 
hull, 


Z. INDX4 is the Index of the point below INDX1 (lower). For the same reason 
as in “1, above, this point is located either | index higher or lower than INDX1 de- 
pending on whether the outer or inner hull is being considered. The resulting expres- 


sion for the index becomes 
INDX4 = INDX1 +1) * (L61) = 1 * (Lé2). 


FIGURE 3a, FIGURE 3b. 


All ship section points lying between OP] and OP2 are transferred to 
the Compartment Section Offset Table. 


K = sentinel Indicating whether OP1 lies on a ship section offset (whose 
index is NWP) or between two ship section offsets (whose indices are 
NWP and NWP+1). The values of K for these two conditions are 2 
and 3 respectively. 


NWP = index of ship section offset on which lies OP1, or just below OPI if 
OP] is between 2 offsets. 


Since only offsets lying between OP! and OP2 are to be tronsferred, the 
index of the upper offset to i transterred Is : 


LMT = NWP - 3 +K. oe 


awe eee 


SUBROUTINE MODIFY (continued) 


For example, in Figure 3a., NWP = 6 and K = 3; thus | 
LMT =6-3+3=6. | 


In Figure 3b., NWP = 6 and K = 2; thus 
LMT=6-3+2=5. 


INT = index of lowest offset to be transferred. LWP is NWP fer the lower 
outline point, OP2. Thus 
INT = LWP + 1. 


In Figure 3a., LWP = 2; thus 
INT=2+1=3., 


In Figure 3b., LWP = 3; thus 
INT=3+1=4, 


In summary, offsets 3 to 6 are transferred in Figure 3a., and offsets 4 to 5 
are transferred in Figure 3b. 


4, INDXS is the index of the outline point under consideration. The program 
encounters this statement when considering either the lowest outline point, INDX2, 
or any other outline point, INDX1. Thus, INDX5 equals either INDX2 or INDX1 
depending on the situation. Since L5 = 1 when the lowest outline point is under con- 
sideration, and L5 = 2, when any other outline point is under consideration, the se- 
sulting expression for INDX5 is 

INDX5 = INDX2 * (2=L5) + INDX1 * (L5=1) 


5. 


Bulkhead A Bulkhead F 


Sots 


Lt | | 
OP, 


\ 


FIGURE Sa. 
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SUBROUTINE MODIFY (continued) 


CO Poe At AN a A, 8 Gm ET Se ARTO! ONE 


Bullhead a ~~) bulkhead F 
Intermediate lee 


ee coe OP ¥ 
OP, longltudinal bulkhead 
FIOURE 5b, 
P bulkhead F 
~Bulkhead A 


~- Intermediate Station 


FIQURE Se, 


In the case where a compartinent has u boundury of a sloping lonyl tudinal 
hulkhead a slopligdeck, the Intermediate station must? be modi fled that a Ilnea 
relationship exlets between It and the bounding compartment bulkheads, 

Thue It (Zop,° Zon f ), the corresponding Z at the Intermediate station 
(Zep Zon)’ arhe)4 ' Zope. 
if (Yop, . Yor) fp 0, the corresponding Y at the Intermediate station ls 


ors” Tope al * Yon, 


SUBROUTINE MODIFY (continued) 


AK 


INDX2 


INDX2¢1 
FIGURE 6a, FIGURE 6b, 
INOX2 ta the Index of the fret lowest outline point encountered, The Inner 


hull normally startet INDX2, as in Pigure 6a, However Ifa» 7 between INDX2 and 
INOX2+1 equals sero, avin Ngure 6b, the Inner hull starte at INDX2+1, 


SUMOUTINE MODIFY 


NPT-NPIS(JSTAT) 


_ 
DO 1000 J-1,NPi ) 


2ST (1,J)-21S(JSTAT, J) 
VIST (1 ,J)-YIS(JSTAT, J) 
VIS(JSTAT, J) - 0 
L\S(JSTAT,J) - 0 


JOPST(1,J)-JUPIS(JSTAT, J) 
JOPIS(JSTAT,J) = 0 


ELOWCHART 


Control characte: © 


NPT - number of points on the ship 
tection, outer hull. 


Transfer 2 and Y points of ship section 
outer hull to temporary storage arrays 
ZIST and YIST, ovter hull, and clear 
ZiS and YiS arrays, outer hull, to sero, 


Store "number of puints" on outer ship 
hull liv NPIST(1), 


Transfer breakpoint locations on ship 
section outer hull to temporary storage 
array JBPST, outer hull, and clear 
array JBPIS, outer hull, to sero. 


{ 


_. *) 


SUBROUTINE MODIFY (continued) 


INDX2+1 
FIGURE 6a, FIGURE 6b, 
INDX2 le the Index of the first lowest outline polnt encountered, The Inner 


hull normally starts at INDX2, asin Figure 6a, However Ifa Z between INDX2 and 
INDX2+1 equals zero, asin Ngure 6b, the Inner hull starts at INDX2+1, 


SUBROUTINE MODIFY 


NPT -NPIS(JSTAT) 


DO 1000 J- 1, NPT 


ZIST(V ,J)-ZIS(JSTAT, J) | 


VIST (1 ,J)-YIS(JSTAT, J) 


YIS(STAT, J) ~ 0 
ZISUSTAT,) - 0 


"1000 


~ (Gontinve ) 


JOPST() , J}-JBPIS(JSTAT, J) 
JOPIS(JSTAT,J) = 0 


FLOWCHART 


Control character - 0 


NPT - number of pointe on the ship 
section, outer hull, 


Transfer Zand Y pointe of ship section 
outer hull to temporary storage arrays 
ZS and YIST, outer hull, and clear 


ZiS and YIS$ arrays, outer hull, to sero, 


Store "number of pointe" on outer thip 
hull In NPIST(1), 


Transfer breakpoint locations on thip 
section outer hull to temporary storage 
array JOPST, outer hull, and clear 
array JOPIS, outer hull, to aero, 


00 10(1030, 1020), ICAI 
2 


41920 
INOX! - JSTAT ¢ LSTAT 


NPIST(2) - NPT 


DO 1030 Je, NPT 


Z181(2, J)~Z1S(INDX1, J) 
YIST(2, J)e¥IS(INDX 1, J) 


VIS(INDX!,J) - 0 
ZIS(INOX!,J) -_ 0 


DO "1040 J-1, 10 


JOPST(2, J)~JBPIS(INDX!, J) 
JOPIS(INDX1,J) = 0 


1) 040 


ICAT - 2, ship has 2 hulls GO TO “1020 
to transfer Inner hull to temporary sturage. 


Index of Inner hull - Index of outer 
hull + number of stations In compart 
ment hull. 


Number of polate In inner hull, 


Store number of polnty on Inner hull 
In NPIST(2), 


Tranefer 2 and Y polnte of statlon, Inner | 
hull to temporary storage arrays 

ZIST and YIST, Inner hull, and 

clear ZiS and YIS, arrays, Inner 

hull to sero 


Transfer breakpoint locations on Inner 
hull to temporary storage array JBPST, 
Inner hull, and clear array JOPIS, 
Inner hull, te sero, 


SUBROUTINE MODIFY FLOWCHART (continued) 


Initiall sation 
(See Symbol List) 


NPT- + NINT(LBULA) 
SQUNWHALE - 9000,0 
XKEEL * 9000.0 
SHELL ~ 9000.0 


IF(ABS(ZINT(LAULK , NPT) 
ZINT(LBULK ,NPT=1))».01)4,4,5 Drop conaideration of laut 
: outline point If Z of that 
polnt ls same as Z of second 


to last outline point. 


Locate lowest 
outline point, 


If lowest point ln at keel, 
lowest paint leat Index J) If 
above heel, It lsat J-!, 
Store Index of lowest polnt 
In INOX2, 


SUBROUTINE MODIFY 


La 
INT 
INTI 
Lé 
Ld 
L4 


NPT- 
XKEEL 


INOX | 
INOX) 


SHELL 


IP(ABS(ZINT(LQULK , NPT) 
ZINT(LBULK , NPT 1)» .01)4,4,5 


t t ' t 4 


| 
2 
2 
2 
we 
- | 
=| 


+ NINT(LBULK) 
GUNWHALE ~ 9000.0 
- 9000.0 


* 9000.0 


IP(ZINT(LBULK , J), OT .XKEEL) GO 7O 9112 


_INDX2 + Jo) 


he 


FLOWCHART (continued) 


Initlaltsation 
(See Symbol let) 


Drop consideration of last 
outline point If 2 of that 
polnt ls same as Z of second 
to last outline polit, 


Locate lowest 
outline point. 


If lowest point ls at keel, 
lowest point is at Index J) If 
above keel, It ls at J-l, 
Store Inden of lowest point 
In INDX2, 


SUBROUTINE MODIFY 


191): 
LMT - INDX 
LMI - INDX2 


VF(VINT(LBULK J). LT, SHELL) 
QO 10 15 


no 


IF(ZINT(LBULK, Jv 1),LE, 
__LINT(LBULK ,J)) GO 10 16 


ne 


FLOWCHART (continued 


LMT - number of aur!ine polets 
for thle hull (Inner or 
outer) of compartment, 

LMT1 ~ Do luop limit. 


IF ICAT = 2, ship has double hull, 
HICAT - 1, ahip has single hull, 


If frat outline point lean shell, 
It must be an outer hull since 
points are taken clockwise, 
OO TO 16) !f not, teat further. 


Tat all outline polite an thle 
compartment hull, 


If point does nct Ile on shell, 
go to teat nent point, 


Goln clockwlas, and assuming 
that once thy bulkhead Intersects 
the shell [t must follow the shell 
for some distance, then the nen! 
outline polni must also be o shell 
Intersection, Then lf Z of the 
2nd polnt ly less than Z of the let 


the Intersection must lle on the outer 


ship hull. If Z of the 2nd paint 

ls greater than Z of the ler, It 
must Ile on the Inner ship hull and 
the ship hull Index, INOX9 must 
be reset to 2, 


iia: 


SUBROUTINE MODITY PLOWCHART (continued) 
6) 
D 
Ne 
ae NPIST(NDX3) | Ninberiotpelntcan wiphull 


i INeINTHLAULK, 44NDX1)¢1.QUNWHALE) Sy 
Pd 


GO 10 18 
no If bulkhead goes io top of thip 
a aane oe ee section, Z value of upper out- 
ZINT(LBULK, K INDX1)- “LIST(iNDX3, NPT 1) line point te sat equal to Z of 
a Paik NS ak ig ay upper polat on ship section, 


Index of cutllne point where 
whetltutlon was made, 


ANB. sowie 
[ saris(istar, 1)-1] Firat paint on compartment 
ase aes tation ls always a brenkpolnt, 


Abe ES me ae 


we GUTO 20. 
VE(ZENTILOULK ,INODX2),O1 XKEEL.) 2 
; no If bulkhead goes to bottom of 


thip section, Zane ¥ of firat 


| ZIS(JSTAT 1)-21ST(INDXS, 1) polnt on compartment seetlon 
YiS(JSTAT,1)-YISTINDX3, 1) | are te! equal to lowest! polnt 
| on thip sectlon, 


Sentinel Indicating first compart 
ment point Iles ona ship point 
(K-2) and Index of that point 
(NWP+] )e 


Sentinel Indicating use of loweet 
compartment gle point. 
(Sen Statement /125) 


6? 


ia oR Sa ae Sang RL aia ee ae Ian AERC Raa 


SUBROUTINE MODIFY FLOWCHART (sontinued) 


Incroae Index of compen tment 
rection polnts by | av storage 
of pulot wae jut completed, 


Sentinel Indloaling use of com 
purtment outline polit under 
vonideratlon (INDEX-!NOX1, 
see wien around Statement #125) 


AMT ) Do lor euch outline polnt on 
thts compartment hull, 


INDXI-(LMT =H 1)* (L6=4) 01 #(26L8) | Index of oviline polny 
a ae ‘ 


(See notes, 


Lb | )=1*(L6=2) Index of nent highest outline 
polnt. 


(See notur, 12) 


INDX4-INDX1+ 1° 


2 RI TI «68 noo 


1F(, NOT, If the current outline point to 


(YINT(LOULK ,INOX1),G2@, SHELL. AND, on the shell, and elther the pre~ 
(YINT(LBULK ,JNOX4),G2, SHELL.OR, ZINT vious (lower Z) outline pulnt 
(LOULK ,JNOX4), LE, XKEEL})) leon the shell or le the heel, 


GO 10 1435 those pointe on the thip section 
lyliy between the two outline 

polnte must be transferred to the 
compartment wectlon; goon, If 
thle statement ls nol true, the 

bulkhead dows nol Intersect the 
shell between these two polnte, 


GO 10 Statement #125, 


SUBROUTINE MODIFY 


[M- NPIST(INDX3) 
CML ay 
[ restaczinanes out, INOX!), 2157) 


Neen) sy snr lS 


[ umt- NWP 1K 


ores OD Bw 


INT-LWP el 


- = + —-0(H0 16 TINT IMT) 


IS(JSTAT, 1.2)  ZIST(JSTAT, J) 
cv Eee Y= YIST(ISTAT, J | 


FLOWCHART (sont!nued) 


Store Indax, NWP, and location sen 
tinel (Indluating polnt lyon NWP or 
betwaen NWP and NWP +1, Gee Sub- 
routine TESTI), K, of last outline polnt 
rhut Hes on shell, 


Number uf pointe on thip hull wetlon= 
for Subroutine T2ST), 


Find talative location of outline point 
In sblp section table of offerte, 


Indax of highest point on ship section 
to be hansforred lo compartment rection 
(See nutes, 99) 


Indox of lowest polnt on ship section 
to be transferred to compartment section, 


Transfer ship section offsets, inden J, 
to compartment section offest table, 
Index 42, 


Increase Index of compartment seatlon 
offeet table, 


om ee, A 


SUMOUTINE MODIFY 


TROOPSTIINDXS,N)=J) \ 
en LPL 


cue 44)-b2-1] 


oD 


\\ 
LJ-Loel 


’ oa 
(1) Pentlnus) 
4) 


‘ 
(Cant nue | 


ZIS(JSTAT,L2)- 21 NT(LBULK ,INOX4) 
+(ZINT(LBULK ,JNDX 1)» ZiNT(LBULK, 
2 


VIS(JSTAT ,L2)-YINT2(ZtS(JSTAT ,L2)) 


ee, |. ei eo, ERS a AE + =: 


FLOWCHART (continued) 


Search chip wotlan break pulnt 
tublo, Many bradakpolnt Index ~ 
Index of affeet ust transferred, 
define compartment sectlon braak« 
point ae inden of uffeet in coms 
partment wetlon atfeet table (4419), 
increate Index of compartment 
tection breakpoint table by | 
(La~Lart), 


Inereave Index uf ship breakpuln 
lable where search will begin, 


If the Incom of the fret hip 
sectlon offiet to be transferred 
In the above DO 120 LOOP Ip 
greater than the Index of the 
last one, no polnts were trant= 
ferred as none lay between the 
two outline points, Therefore, 
an Intermediate offset must be 
created In order to preserve the 
curvature of the sectlon, It ls 
tuken at a point halfway (2 
direction) between the two out 
line pointe, and stored In the 
compartment tection offset table, 
The Index of thle table (L2) Is 
then increased by |, 


eo mew + 


SUBROUTINE MODII'Y 


| INDX INDAY (de UM)OIND (a i: 


IFUMOD FQ 1,.OR,(NINE.NELOAND,. 
(INDXS 002, 1,08, INDXS, EQ, NPT))) 


7 \N(ABS(ZINT(LBULK « 
ena IMOX3)}=..1) 


came 8 
Shad ome \ 
‘ 
ee a 
eae eines 


Oh) Wome see 
CG) us Ustar, 12)-ZINTLAULK, INDX: 


or) oe. 5): } omens 


AIS(ISTAT 1.2) (ZINT(LBULK +) ENIX eo 
ZINT(L BULK ,LNDX5))/(XBLK (I BULK +1) 
= KBLK(LOULK))*(XI S(JSTAT XBL K (L BULK)) 
VZINT(LAUL kK, )|NDXS) 


00 am 6 MAE Seem ane =n eee 


HLOWCHART (continued) 


Index of autline polit, 
See notes, 94, 


Vf thie le @ bulkhead (UMOD~1), 
ut If an intermediute stetion and 
the wutline pulnt in question lo 
the top of either the Inner or 
cuter comparinient hull end alto 
live un the top of the ship seatlon, 
then 2 of the compartment offeet 
ls the 2 of Hp outline polnt. 
(Statement 9136), Otherwise,.., 
(for Intermediate stations unly), 
ont A Z Letwoon the outline 
pulnt an the hulbheads fore and 
alt of the Intermediate statlun, 

lf A 290, the Z of the compart~ 
ment offee? lo the Z of the vutline 


polnt. 
(Statement 4136) 


\f £424 0, Interpulate the Z 

of the compartment offset by a 
linear relationship between the 
Jiaf the outline points on the 
turroundling bullhead, 

(Sue nutes, 95) 


«0 


SUBROUTINE MODIFY 


| IMVINTELBULK I NDKB).LT. SHELL) S, 
: 09. JO 14. 


Ve | 
VIS(STAT,L2)-YINT2(ZIS(USTAT 42 


OO TO (29, 490),L5 


’ 
GO TO (175,180), JMOD 
# 
Yi S(J STAT L2)-YINT(LBULK ,|NOXS 
/ 
VP(ABS(YINT(LBULK+ 1,1 NOXS)- 
YINT(LBULK ,ENOXS)> ,1)190, 190,200 


id 


VIS(ISTAT, b2)-ViNT(LBULK,INDX3)] 


VIS(ISTAT ,L2)>(YINT(LBULK+! ,INOX3)= 
YINT eOULK ,JNOXS))/ (XOLK(LOULK + 1 )= 
XBLK(LBULK))*(X1S(JSTAT)=XBLK(LBULK)) 
+YINT(LBULK ,INOXS) 


PLOWCHART (continued) 


It the outline point te within 
the shall, OO TO “145, 


If 1) Vee on the shell, the ¥ 
value af the compartment offset 
- the ¥ value af the chell a 
the Z of the autline point (vee 
Interpolation routine YINT2), 


IP L3D-1, outline polnt te lowes 
point, OO 10 439, If Ld-2 
In le not, GO TO 450, 


if thie le @ bulkhead (JMOD-1), 
the Y value of the compartment! 
offvet - the Y value of the out= 
line polnt, 


(This le an Intermediate station. 
and the outline point! Iles within 
the shell), Tes A Y between 
the outline point on the bulkheads 
lore and aft of the Intermediate 
station. If A Y«O0, the Y¥ of the 
compartment offset ls the Y of the 
outline polnt (Statement 190), 


I! AY 4 0, Interpolate the Y 
of the compartment offeet by a 
linear relationship between the 
Yo ol the outline polnts on the 
turrounding bulkheads, 

(See notes, 45) 


SUBROUTINE MODIFY 


Ao. Ub: end aoe o8 
nis S(JSTAT Lapeb 2 


i 


[ua-440| 


a 


FLOWCHART {eontinved) 


Sore index ul compatinent 
offee' In compartinent break 
pulnt table, (offeet |uet stored 
waa outline point and all auto 
VVever pwolevtn erie boca gaalinte), 


Ineroare compartment! break pint 
lable Inedon iy |, 


Iner@ase cuinpantinent « Hee 


lable Indew hy |, 


If the outline pulnt jus’ opernied 
on lay within the shell, go on 
to the nent outline point, 


. VP(VINT(LBULK,INDX1), LU. HALL) > 
9010900 


ye 


[cau TPSTO(ZINT(LAULK,INDKI) 2151) It Ht lay on the shell, find the 


relative location of the aulline 
polnt In the thip section table 
of offsets, (M le the number of 
pulrio un the whip ewctlan). 


~ 


Pao RAI I BTS OO ISADORE 


| ; — 1g Gontinun Go an to nent vutline palit, 


[aa crnaiy 


-_ Oe ee cre. 
< (NINE, Nf, O)ZINT(LBULK, Nil ) 999, 2 


The number of paotnte on the 
compartment sation, 


NINE 4 0, the top ostlin 
polit liv the wutlon lay on the 
tup of the tentlon and the autiuul 
Z vale replaced the slynal 
value of 9999, 0LSee Staten 
"1601, pagaD ). Restore slynal 


; Rf value of 9999.0, 
wo) 


SUBROUTINE MODIFY 


340), JCAT 


Aah a NERD § Sar -O 


OO TO (863,510, 


f 
€ GO 10(580,510,),ICAT 
ech TONNE IE H y oiuanicicaes 


JSTAT-JSTATILSTAT 


|NOX!-NPT 


IP(ABS(ZINT(LBULK , INDX2)- 
ZINT(LBULK ,I NOX2¢1))- .01) 
330,390,520 


INDX2-INDX2¢1 


LLOWCHART (continued) 


1 ICAT - 1, bulkhead outline polit 
definittune Indloate alngle hull; QO 

TO "505 to 100 te hip hee double hull. 

IE JCAT - 2, compartment has double 
hull, outer hull here [vet been mod! fled, 
GO 10 9510 tor Inner hull Initlalleation, 
If JCAT - 0, compartment has double 
hull, and Inner hull hee fuet been madi « 
Ned, OO 10 4340, 


IICAT - 1, ship hee eingle hull, OO 
10 "990, IHICAT - 2 ahip hae double 
hull and therefore compartment must aleo 


have double hull; GO TO 4510, 
Signals modifleution of Inner hull, 


Index of Inner hull - index of outer hull 
‘number of compartment stations, 


Index of outline polit defining top of 
Inne: compartment tection, 


It 2d Z - 0 between lowes! outline point 
and next outline point, Increase Index 


‘ of lowest outline polnt by | (to second 


point) (See Notes 6), 


70 


Pee eS ee a em oe eee, eee eee, Si ie Ea  .sbiia.s.. Ab © > 


SUBROUTINE MODIFY 
¢ 


INTY = INDXa+ | 


INT = INDX2 
LMT ~ NPI 
LMT)- NPI 
lo + | 

La | 

ld - 2? 
4-2 


Pree oo eo 8 


029 1564 °K" ) 
[NPT-NPIS(STAT)] 
DO 543 1-1,NPI ) 


INDX1t-NPT=141 


IF(ABS(VIS(JSIAT,INDX1))=.01) 
TORT LL, ee 


Continue 


NPISJSTAT) ~ | 


FLOWCHARS (continued) 


Inltlaliaation far inner compar t= 
ment hull, 
(See Symbol 111) 


HO VCAT ~ 1, oblge hee olngle hull, 
OO 10 16, IFICAT - 7, ship 
hen davble hull; OO TO 14 to oat 
ship hull Index to Inver ship hull 


Kew! JCAT, 


Number uf points on Inner come 
partment hull, 


Search all pulnte on Inner hull, 


Start at highest polnt and work 
down (work backwards from end 
vf tabla) 


Ht the Y af the offvet ls aero look 
at the nent point) If YIS 4 go to R, 


If all offsets have sero Y, there ls 
effectively no Inner hull, Set 4 of 
points on section Inner hull to |, 


-m -RBOGs. & tm 2. = 


— — 


SUBROUTINE MODIFY 


TRU.EG, 1) 990,546 
<i ae ae. 


IM(ABS(ZISUUSTAT,INDX1) 
= Z1S(JSTAT, INDX 1+ 1))> .01) 
pee Lae LI6CL, emer 


25 OBA Ee OR wee en HE a | OUR 6 RL IPN 


Esa INOX? | 


| 
tes be oan-0) 10 0h PANS CU Ne 


[ NP S(JSTAT)“INOX 1 #1 


Ema t ES = cee 


FLOWCHART (continued) 


If offeet haa a Y of sere (I-1), 
yo to “950, If there lt at least 
one offret witha Y of sera go 
to "946, 


If > 2-0 between the firet point 
of positive Y, INDX1, end the 

laet polnt of sero Y, INDX1+1, 

the leet valid point en the stetlon 
leat INDX1, end therefore, the 
number of points on the Inner 
compartment section equals INOX1, 


( ZXZ-0) = the lav valid point 
leat INDX 141, 


Set all Z above last valid poln 
to aero, 


Return to Subroutine GENERATE, 


7\ 


SUBROUTINE DAMAGE 


1, PROGRAM DESCRIPTION: 


A. General 


This subroutine calculates the damage stability characterletics of the ship minus 
selected flooded compartments for any draft, ial and trim candition, There 

Is no balance of the ship); the properties are calculated at th set values of draft, 
heel and trim, The compartments aro subenittad In the same manner as the ship, 
namely, by a series of offsets ona minimum of J stations. All compariment must 
be submitted and thelr ansoclatad Compartment Date Table stored on the Ship 
Data Tape (LU 30) prior to execution of thls routine, 


B, Basic Logie 


1. Input Seation: 
a, The number of compartment sete (NSET) le read In, 


b, For each set of nompartments the number of compartments In the set 
(L), thelr eorlal numbers (L2COMP(!,J)) and permeabilities (PERMB 
(1, J)) are raad In, 


o, The number of Input trime (N4TRIM), drafts (N4DRFT), end heels 
(N4HEEL) are read In, 


1ON4TRIMY 0, then the fret trim = DITRIM()) © Ie set equal to sero, 
Otherwlee, the Input trims are read In and the frst ene-DITRIM(1) 
Is wet equal to sere, 


a 


The Input dratte (DIORFT(1)) are read In, 


<= 
- 


The Input heel angles (DIHEEL(I)) are read in, 
g: A llatof the Input data ls printed out, 
2, Caleulation Section: 
a, Tha Inter-program sentinels KW, ITEST!, and KW ore defined, 


b, JSERNO Is defined as 
ISERNO, the serlal number of the ship, 


6, The set Index | le Initiallaed, 


d. The trim Index ITRIM Is Initlallaed and the trim angle converted 
from degrees to radians, 


@, The draft Index IDRFT Is Initlallaed. 


The hwel Index 13 | Initlaliaed and the heel engle converted from 
degrees to randlans, 


—_ 
- 


ISERNO, the sarlal number of the ship or compartment belng oper= 
ated on, le set equel ta JSERNO, the serial number of the ship. 


The Subroutine OFFSET ls cullad, where the Ship Dete Table ls reed 
Irom the Ship Date Tapa, 


= 


The Subroutine VOLUME ts called where the volumetric proper!ies 
of the ship ore calculated, 


The volumetric propestio: of the ship are stored In VISHIP, TCBISM, 
XLCBISH and VCAI SH, 


zx 


L ls net equal to the rumber of compartments In set |, 


The compartment Inden 14 le dofined, 


1SERNO |e set equal ta the serial number of compartment 14 Inset |, 


3 


ITEST! Ie oe) equal to 2 (this signols OFFSET that the Compartment 
Data Table ls on the Ship Data Tape), 


_ 
—_ 
- 


(OFFSET Is called whare the Compartment Data Table le read from 
the Ship Datu Tape. 


p. VOLUME Ie cal vd where the volumetric properties of the compartment 
| are calculated, 


The volumetric propertion of the sompartment are stored In VOLIC(I4), 
TOMIC (14), XLCAIC(14), and VERIC(I4), 

If all compartments in set | have been calculated, go to S) otherwise 
Inaremen! comparimant Index 14 and go tom, 


2 


= 


6. The compartment tial aeeumulators are set to vero, 


1, The volumetric propertion uf the compartments in tot | are summed 
and thelr Individual values printed Gut, 


The volumetric properties of the tntal flooded volume and the Intact 
thip are printed ou!, 


c 


The flooded volume |s subtracted from the Intact ship and the net 
properties of the flooded ship are printed eut, The righting erm le 
computed and printed out, 


< 


w. Inerement heel Index 19 and yo back tog) Ifall heel angles have 
been calculated, Increment draft index IDKFT and go back to f) 
If all drafts have been calculated, Increment trim Index TRIM and 
go back toe; Ifall trims have been caleuvlated, Increment compart= 
ment te! Index | and ga back to d) If all sets have been calculated, 
return to executive program, 


INPUT 
A, COMPARTMENT DESCRIPTIONS 
A.! Method |= Submit Compartment Offrets 


Zach compartment Is submitted ava regular "sip" aacording to tho Instructions 
pertaining thereto, For each compartment a minimum of 2 statlans are required. 
They should be submitted In accordance wth date cards 1, 2, 9, 4, and §, 


Card Type 'i CC 48-30 should he punched accordingly: 000) com= 
partment on both alder of centerline) b+! comeai:~ 
ment on starbourd side; b= !, compartment on port alde, 
(b represents blank), 


Two restr'etions enlat at present on the method of compartment deserlption which 


muel wi lowed! 


1, The station spacing used In describing the compartments must be Identical to 
tha? of the parent ship, |.e., the station numbers on the Type 3 Data Card must be 
In terme of the parent ship, 


2, The "LBP" on the Type 3 Dota Card must be that of the parent thip, not the 
length of the compartment, 


When the compartment description data lt submitted the Type 2 Data Card should 
be blank © thie signified that no actual major wpervitors are tn be executed for the 
compartment) only the Data Table will be formed and stored un tape, 

A.2 Method 2= Submit Bulkhead Ourline = See Subroutine CMPRTMNT 


The compartments are Identified to the damaye stability routine by the serial 
number that appeared on the Type | Data Card used in thelr submission to the program, 


Once the compartments are submitted they are stored on the Ship Data Tape and 
need not be resubmitted for future runs, 


B. INPUT FOR DAMAGE STABILITY ROUTINE = See System Description 


Il Output 


Table of volumes of fluoded compartments VOLIC, 
Vertleal centers of buoyanay of Hooded compartments = VCBIC, 
Transverse centers of buoyancy of flondad compartments = TCBIC, 
Longitudinal centers uf buayuncy of fusded compartments = XLCRIC, 


Total flouded volume = VOLIT, 

Vertical center of tuoyancy of trtal feuded volume = VCRIT, 
Traneverte center of bunyancy of total flncded volume = TCBIT, 
Longltudinal center of husyane, of total flooded volume « XLCBIT, 


Volume of intuet ship = VISHIP, 

Vertleal venter uf buoyaney of Intact ship = VCBISH, 
Traneverte center of buoyancy uf Intact ship = TCBISH, 
Longitudinal center of buoydney uf Intact ship = XLCASH, 


Volume of damaged ship = VOIIT, 

Displacement af damaged ahin « SPL, 

Vertical center of buoyuney of damayed ahlp = VCBIT, 
Traneverse center of bunyanay af damaged thip = [CHlT, 
Long!tudinal canter of buoyunry of dumuged ship = XLCBIT, 
Righting arm of damayed ship » ROTAKM, 


NOTE: Finaded volumes are ustual volumes « Input permeabilities, 
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FUNCTION YINT2(TEST) 


This program Is identical to WPLANI except for the following Items: 


1) YINT2(TEST) ls called as a function subprogram rather then as a sub= 
routine, 


2) The tes: value of draft at which the halfbreudth ls desired le now Includ= 
ed In tha calling statement as TEST, rather than stored In COMMON, 


3) The following array ond varlables have been replaced as Indicated* 
(for deser!ption of the meuning of the different arrays wee the documentation | 
for subroutine MODIFY): 


Zi ST---pZ | 
YIST---pY 

YINT+-<@ YI(ISTAT) 

Le==-pM 

\BREA=-~p) BREAK 

INOXIJ-<=pi STAT 

JOPST+=--=p JREAK 


4) The function FINTP2 ls now written out explicitly. 


5) Por further detalls sew NavShips 0900-006-5810. 


*The symbol "==" means "replaces". 


SUBROUTINE TESTI 


This oubroutine Is slmilar In loylc and objective to Subroutine TEST! of the 
original Shia Hull Characterletic Program, namely, the determination of the rel= 
ative postion of a test value with respect toa ‘able of sequential numbers. It le 
uied for locating bulkhead outline polnts on the ship section, 


Two sentinels are OUTPUT: 


Ki Aentinel indicating the status of the teat value with respect to the array 
defined follaws: 


K ~ | the value of TEST le lees than the maximum value In the array, 

K «2 the value of TEST = a value In the array, 

K ~ J the value of TEST Iles between 2 values in the array (Indloes 
NWP and NWP+1), 


NWP: The Index In the array of the value In array corresponding to TEST(K-2), 
or the first value less than TEST(K-3), 


Pur further detalls see Nav Ships 0900-006=5750 where the following arrays 
and varlables have been replaced as Indicated: 


NPT=-=»M 
INOX3=--@/STAT 


PROGRAM MODIFICATION TO SHCP (YGHULL) 


| Por Double Offset Modification. 


A, The following changes must be made to DIMENSION variables In all 
programe: 


All varlables with @ subscript 1S must have the flret dimension which 
18.30 changed to 60. 


All varlables ending In 28 must have their dimension changed from 
30 to 60, 


Ld 


The following draft variables must alto have thelr dimensions Increased 
from 30 to 60: H, HIHOG, HIMAR and HISAO, 


&D 


The flrat dimension of JOREAK shall be Increased fram 30 to 60 and the 
flrot dimension of the variables $ shall be Increased from 30 to 60, 


> 


The flrat dimension of the varlables Y and 2 shall be Increased fram 
IC to 60, 


a 


the dimension of the varlables X, NP, STATNO and YI thall be In | 
creased from 30 to 60, 


o 


7, YIINT(4,30) shall be added to DIMENSION, 
B, The following variables shall be added to COMMON: ICAT, LIMIT and 
YIINT, 


The following lea Ilet of required program modifications, referenced by the 
card numbers that oceur In the program Ilatings In the Individual program 
write=upe, 


C. Hull Characteristics Bxecutive Program Nav Shine 0900-006-5600 


INVSYM,ICAT 
ILO, XMID, JOUMMY , LIMIT 
1920 DO 1141 1STAT#) LIMIT 
ILeP, XMID , JOUMMY , LIMIT 
DO $462 (STAT+1,LIMIT 
IP(ICAT=2)728, 1289, 1269 
728 KW-! 
1289 IPROG=! 
NWISYM,ICAT 


D, Subroutine OFFSET Nav Ships 0900-006-5720 


DC 1756 ISTAT«1, 60 
16P,XMID, JOUMMY LIMIT 
JI9 DG 390 ISTAT$1 LIMIT 
140 OC) 141 ISTAT~1, LIMIT 
ICAT-JOUMMY 
VP(ICAT+2)999 994,395 
334 ICAT#J 
ISTAT-ISTAT#1 
LIMIT~2*NSTAT 
60 TO 10 
335 ICAT~2 
NSTAT-NSTAT/2 
GO TO 346 
339 LIMIT+NSTAT 
346 DO 411 ISTAT#1 LIMIT 
OO 65) ISTAT=1,LIMIT 
1L OP, XMIO ,ICAT LIMIT 
DO 190 ISTAT-1,LIMIT 
DO 660 ISTAT~1,LIMIT 


E, Subroutine HYDRO Nav Ships 0900-006-5700 
ILIMIT-LIMIT 

890 LIMIT-ILIMT 
RETURN 

F, Subroutine AREAS Nav Ships 0900-006=5700 
1a-1 STAT 


GO 1O(3110, 889,690) ,|CAT 
809 VIINT(1,1STAT)*“¥98(1) 
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HullOalo 
Huli0680 
Hull0870 
Hull0940 
Hull0950 
Hull 1610 
Hull 161) 
Hull 1860 
Hull3990 


OFFS0740 
OFFS0830 
OFFSO910 
OFPS0940 
OFFS$099) 
OFFS177| 
OFFSI772 
OFFSI779 
OFPPSI774 
OFPS1775 
OFFSI776 
OFFS$1777 
OFFS1776 
OFFS1779 
OFFS1780 
OFFS2110 
OF F$2400 
OFFS2410 
OFFS2690 


bYOR 1681 
HYOR1860 
HYOR1061 


AREAO799 
AREAII00 
AREAIIO! 


00 70 JI10 AREA3I02 


890 I1-ISTAT-NSTAT AREAII09 
YIINT(9,11)-¥58(1) AREAII04 
GO 10 9110 AREAIIOS 
VF(ICAT=2)777 778,779 AR2AI964 

776 $\-S2S(ISTAT) AREAIIGS 
AREA-AREAQ $(/ STAT) AREA3386 
1CG-TCG28(I STAT) AREAII07 
VCG-VCG2S(I STAT) ARIAII68 
| P(ABSP(THET1H)= 0001 )667 , 667 , 668 AREAIIB9 

668 YIINT(2,1STAT)-Y58(2) AREAI290 

667 ICAT~J AREAII9) 
ISTAT-ISTAT¢NSTAT AREAII92 
GO 10 222 AREAII99 

779 VEGIS(I2)-(VCG*ARBA=VC G2 S(ISTAT)*AREAQS 

(ISTAT))/ =AREAQ S(1S AREAII94 

\TAT)) AREAII9S 
1CG2S(12)-(TCG*AREA=1CG28(i STAT*AREAQS(ISTAT))/ 

(AREA+AREA2S(IS AREAIIV6 

\TAT)) AR2A9997 

ARCA2S(12)-AREA-AREA28(/ STAT) aREAIIVG 

$2S(12)-S1+S2S(ISTAT)=2 0° S(ISTAT,L) AREAIIVP 

ICAT~2 AREAI400 

I P(ABSF(THET 1H)» 0001 )664 , 664 , 665 AREAI401 

665 YIINT(4,1STAT)~¥38(2) AREAI402 


664 STATI STAT-NSTAT AREAI403 


PROGRAM MODIFICATION TO SHCP (YOHULL) 
(continued) 
\! Modi fleations for incorporation of Subroutines DAMAGE & CMPRIMNT 


The following |: a Ilst of required modifications to SHCP, referenced by the 
card number that ocour in the program Ilstings In the Individual program write-ups: 


A. Hull Characteristics Gxesut!ve Program Nav Ships 0900-006-5600 


18°, XMID , JOUMMY LIMIT, KOUMMY HULLOBSO 
ILOP, XMID JOUMMY , LIMIT ,KOUMMY HULLO940 
62 GO TO (20,70,90,3999 , 3999, 1142, 40,95) MULL 1890 
95 CALL CMPRIMNT HULL3775 
GO TO 60 HULLI776 
100 GO TO (390,340,350 360,370 , 380,390,500) ,/TESTI HULLI820 
390 CALL DAMAGE HULLI97 1 
IPROG-=IPKOGs1 HULLI972 
GO 10 140 HULLI979 

B, Subroutine OFFSET Nav Ships 0900-006-5720 
18P,XMIO ,JOUMMY , LIMIT, KOUMMY © FFS0600 
1)S¥YM*KOUMMY OFFS0992 
ILOP,XMID,ICAT, LIMIT, LISYM OFPS2400 
WRITE(30) IEOP, NAME, NSTAT ,NP,XLOP,XMID OFFS2470 


ICAT,LIMIT,IISYM 


DATA PREPARATION 


The required Input for the Hull Characterlstlos Program depends upon which major 
whroutines are desired to be executed, The Card Type Usaye Map on the following 
page Indicates the required card types for each of the major subroutines, A partloular 
card type may require more than ono card and thle fact ls Indicated in the Input det= 
cription, 


For any one run of a thip only one of each pertinent card type may be submitted, 
Start at the highest number program being called and work hack, adding these card 
types required for ather called programs not required by any higher called program, 


Only one card Type E ls required per run, regardless of the number af ships belng 
run, It le placed at the end of the data deek, 


Cards type 9, 10, and 14 are optional = see card type desuriptions, 


Data cards for any computer run are submitted in ascending order by card type 


number, axncept for Card Types 25 through Jo, See Systam Dusariptlon, pages 2 and J 
| of this report, for detalls, ie Scare! 
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MACHINE REQUIREMENTS AND CONTROL CARDS 


Thle version of the Ship Hull Characteristics Program Is written In Fortran IV 
for exeeution on the Contrel Date Corporation 3600 System, 


Core required for the system Is 32,590 words (decimal) consleting of 16,686 
words for semmon storage end 15,044 words for the programs, 


To run the system on other than COC equipment, compare the following Fortran 
statements In COC Portran IV to these on the considered equipment, 
1, DATA 


2, COC Fortran IV permits variable name up to @ characters: the only cases 
of names longe: than 6 characters in the system ara Subroutine GENERATE 
and its call in CMPRIMNY, Subroutine CMPRTMNT and Its call In YOHULL, 
and the variable GUNWHALE In Subroutine MODIFY, 


The ecard check set up for running on the COC 3600 |e Illustrated on the follow- 
Ing page. 
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